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Reliable measurements of cosmic dust abundances have been obtained by ionization 
detectors during particle impact and by collectors controlled either by inflight shadowing 
or by penetration-hole identification. A description of the techniques used is given. 

Crater-number densities observed on the lunar surface and on lunar samp es repre 
sent an important source of information on cosmic dust fluxes. The related resu ts from 
the Apollo 11 and 12 missions are reviewed. The overall knowledge gained from these 
measurements leads to the following flux model: The cumulative flux * w mass m 
follows the extrapolation from larger meteoroid-size range ( Watson s Law) and can be 


described by 


& = 10~ 14 ^m” 1 [m~ 2 s *] 


The Pioneer 8 dust experiment and lunar samples indicate a depletion of the flux 
at approximately 10 8 g. However, cosmic dust particles exist in interplanetary space 
at least down to 0,3 p diameter. They are interpreted as nonmetalhc particles m e 

The atmosphere shows an enhancement in particles of about one order of magnitude 
compared to the flux in interplanetary space at 1 AU. No depletion or cutoff could be 
detected. These particles are interpreted as lunar debris or as disintegrated products 


f The numbers of large lunar craters ( >llfim diameter ) in Mare Tranqmllitatis and 
in Oceanus Procellarum are com pared with the meteoroid flux. These comparisons lea 
to a time-variable flux of *-e r», with B = 2.6 and t = time tn 10 9 yr. Thus, the mete- 
oroid flux at the formation of the lunar maria was approximately 4 orders of magnitude 
higher than today. 


T he results of cosmic dust research are 
reviewed in this paper. A description of a 
detector based on ionization during particle 
impact and of a controlled collection technique is 
given. Meteoroid fluxes today and in the earl\ 
solar system are formulated from the distribution 
of craters in various areas of the Moon. Finally, 
the fluxes in interplanetary space, near the Earth, 
and in the upper atmosphere are discussed. 


DETECTORS AND COLLECTORS 

Experience has shown that detectors based on 
ionization during particle impact can achieve 
reliable results. Figure 1 shows a schematic 
diagram of the detector designed for the HE08 A2 
satellite (Hoffmann, 1971). A particle hitting the 
hemispherical gold or tungsten target produces a 
plasma cloud during impact around the impact 
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point. This plasma cloud expands radially. A 
concentrically mounted and negatively biased ion 
collector records the ions from this plasma cloud 
as the electrons are collected at the hemisphere. 
From one impact, the hemispherical target and ion 
collector sense pulses that are identical in pulse 
height, but different in polarity. From these 
coincident pulses one can derive mass and velocity 
of the impacted particle. 

Figure 2 shows schematically the detector for 
the Helios probe (Griin, 1970). The target con- 
sists of inclined strips of gold or tungsten. Elec- 
trically biased grids in front of and behind the 
target separate the electrons and positive ions 
and produce coincident pulses on the grids. The 
ions are accelerated into a time-of-flight tube. At 
the end of the tube the ions are recorded as a 
function of flight time using a particle multiplier. 
Thus, a mass spectrum of the ion pulse is obtained. 
From these measurements information concerning 
chemical composition of the projectiles can be 
determined. 



Impact Spectrometer 

Figure 2. — Helios dust-detector schematic diagram. 


These detectors have been calibrated using a 
2-MV high voltage Van de Graaff electrostatic 
accelerator. Figure 3 shows the calibration of the 
detector in the mass range between 10~ 14 and 
10~ 9 g for iron, aluminum, and carbon projectiles 
in the velocity range between 1 and 40 km/s 
(Hoffmann, 1971). For these ranges of measure- 
ments the following relationship holds: 

( 1 ) 

where (7= charge-pulse height 
yn~ projectile mass 
v= impact velocity 

As shown by Auer and Sitte (1968), Griin (1970), 
and Hoffmann (1971), the rise time of the charge 
pulse — a function of the expansion velocity of the 
plasma depends on the impact velocity of the 
projectile, i.e., on the plasma temperature. In 
figure 3 the rise time of the charge pulse is plotted 
as a function of impact velocity. Further experi- 
ments showed that rise time is independent of 
projectile mass. 

Thus, particle mass m and impact velocity v 
can be determined from the knowledge of pulse 
height and rise time of the ion or electron pulse. 

Ions are produced according to the relationship 
expressed in equation (1) , a function of the impact 
velocity. Not only the intensity of the ion pulse, 
but also the composition of the ion pulse, i.e., 
the mass spectrum, depend strongly on the impact 
velocity as shown by Friichtenicht (1964), Griin 
(1970), and Dietzel (1971). Figure 4 show's an 
example of the spectra for various impact veloci- 
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Figure 3, — Calibration characteristics of dust detectors, (a) Charge pulse per pro- 
jectile mass vs impact velocity, (b) Rise time Ia vs impact velocity v. 


ties of A1 projectiles impacted on an Li-coated 
tungsten target. For low-impact velocities readily 
ionized alkalis can be detected that are present in 
the projectile or arc impurities in the target 
material. Therefore, the Li ions are much more 
abundant than the A1 ions because of different 
ionization potentials. At higher impact velocities 
(5 km/s) A1 ions are as frequent as Li ions. At 
hypervelocity impacts (8 km/s) the A1 line 
exceeds that of the Li line considerably (Dietzel, 
1971). 

The first version of the described detectors has 
been flown several times successfully aboard a 
rocket payload to detect cosmic dust in the upper 
atmosphere (Fechtig, Fcuerstcin, and Rauser, 
1971; Rauser and Fechtig, 1972). Detectors based 
on the same technique are in Earth orbit on 
OGO III and in lunar orbit on Lunar Orbitcr 
35 (Alexander et ah, 1971, 1972), Also, Berg and 
cow T orkers (Berg and Gerloff, 1971; Gerloff. and 


Berg, 1971) have flown a similar detector aboard 
the Pioneer 8 and 9 deep-space probes. For the 
velocity measurements both groups use thin films 
in front of the detectors at the beginning of a 
time-of-flight tube. A first pulse is produced by a 
particle penetrating the film. The second pulse is 
obtained by particle impact on the solid target at 
the end of the time-of-flight tube. Thus, the 
impact velocity is measured. However, there may 
be undesired effects such as particle fragmentation 
(Grim and Rauser, 1969) during penetration or 
instrument cutoffs, especially for low-density 
projectiles at low relative impact- velocities. 

The collectors of the first generation experiments 
(Hemenway and Sobernian, 1962) applied the 
low-angle shadowing technique in the laboratory 
shortly before and after flight to mark con- 
taminants by respective shadows. However, it 
appeared that the contamination problem was 
not solved by this technique. The last mounting 
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Figure 4. — Ion-pulse mass spectra. Projectiles: Al par- 
ticles. Targets: Li-coated Au targets. 


before flight and the first after flight were done 
before the correct shadows could identify the 
contaminants. For example, most round particles 
were produced by the evaporation process itself. 

Recent collection experiments therefore apply 
an in situ shadowing technique, the so-called 
“inflight shadowing” (Skrivanek ct al, 1970; 
Hemenway and Hallgrcn, 1970). By this tech- 
nique the collection surface is shadowed several 
times during flight in different directions. There- 
fore, by different shadows one can discriminate in 
altitude and also exclude pre- and post-flight 
contaminants. Particles that may be produced by 
the inflight shadowing device can be identified 
chemically. 

Another control technique is based on the 


penetration-hole identification first proposed by 
Yaniv and Shafir (19G7), and, in another version, 
applied by our group (Auer et al, 1970) : Particles 
must penetrate a thin nitrocellulose film before 
impacting on a solid target. Figure 5 shows the 
collector (metal plates) with a thin film mounted 
in front. Figure 5 also shovrs a penetration hole 
with a particle-strewn field beneath. The particle 
obviously was a low-density one that frag- 
mentated as it penetrated the foil. 

The inflight shadowing technique has the 
advantage of collecting very small particles 
(0.1-^-diametcr range) at low relative velocities. 
However, the collection area is comparatively 
small (~100 mm 2 ). The penetration-hole iden- 
tification method is applicable for large surfaces 
(~1Q0 cm 2 ), implying a cutoff because of limited 
penetrability of micron-sized particles at low 
relative velocities. 


LUNAR CRATER DISTRIBUTIONS AND 
EARLY METEOROID FLUXES 

The distribution of craters in various lunar 
areas is an important source of information con- 
cerning the distribution of meteoroids in the 
interplanetary space at 1 AU. Figure 6 shows the 
distribution of lunar impact craters as a function 
of crater diameters for the Mare Tranquillitatis 
(MT) and the Occanus Procellarum (OP) as 
published by Shoemaker et al (1970). The large 
craters on the production curve ( > 140 m diam- 
eter in MT and >50m diameter in OP) follow 
a D~ 2 Q distribution law which corresponds well 
with the overall m~ l distribution of the meteoroids 
( D = crater diameter, m = meteoroid mass). 
The deviation from —3 is explained by a variable 
D/d ratio of crater diameters to projectile diam- 
eters. Therefore, these crater numbers arc con- 
sidered as the time-integral of impacts since the 
existence of these lunar areas. Since MT is older 
than OP, MT has more large craters than OP. 
MT averages 2.4 times as many large craters as 
OP (Shoemaker et al, 1970), although MT is 
only 10 percent older than OP according to 
dating results of the related lunar samples (Albee 
et al, 1970). 

A similar relation exists between the number of 
large craters in the Lunar Highlands compared to 
the related number in MT or OP. According to 
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Figure 5— (a) Controlled collector by hole identification, 
(b) Hole-fragmented particle. 

Hartmann (1966} the number of large craters in 
the Southern Lunar Highland is as much as 30 
to 50 times higher than in MT. 

These crater-number densities cannot be inter- 
preted on the basis of a time-constant meteoroid 
influx. Baldwin (1969, 1970, 1971) and Hartmann 
(1970) have plotted the crater-number densities 
of various maria as a function of age normalized by 
the known ages of MT and OP. We tried to 
combine these observations with meteoroid fluxes 
known from direct measurements (Bloch et al., 
1971). 

Let us first assume a general flux function 
<p{m) in differential form or 



#(m)= / (p(m) dm 

J mo 

in integral form known from direct observations. 
If one assumes that the decrease of meteoroids in 
the past is only given by impacts on planets and 
their moons, one can w r rite: 



(2) 


Therefore, the time dependence is 

N=N«r Bl 

an exponential function. For the meteoroid flux 
in the past one can define a flux function /(m, t) 
in differential form or F(m, t) in integral form 
as follows: 


For today t = 0. The flux functions are /(m, 0) = 
<p(m) and F(m, 0)=#(m). To determine B one 
can compare the numbers of craters in MT and 
OP in relation to the related ages: 




<p(m) *e Bt dm dt 


<p(m) *e~ Bt dm dt (4) 

with R = 2.4, the average ratio of large craters in 
MT to OP. From these determinations one finds 

£ = 2.6 (5) 

The time variable flux function therefore reads : 



/(m, t) = <p(m) *e Bt 




F{m } t) = 4>(m) *e 


F(m, t) =<f>(m) *e~ Bt 


(3) 


( 6 ) 
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Figure 6. — Size frequency distribution of lunar craters at 
the Apollo 11 and 12 landing sites according to Shoe- 
maker et al. (1970). 

If one applies this formalism to the model age 
of the Lunar Highlands (LH) of 4.6 -10 9 yr, the 
Lunar Highlands would show 12 times as many 
large craters as MT and therefore the production 
curve would begin at about 3 -km diameter. 
Table 1 shows the enhancement of the meteoroid 
flux at the formation of MT, OP and LH by 4 
to 5 orders of magnitude. 

So far no special flux function $(m) for the 
present meteoroid flux has been used. Gault 
(1970) tried to compare the densities of large 
craters in MT with flux rates based on a time- 
constant influx. He used fluxes according to 
Naumann (1966) and Hawkins (1963) and 


Table 1. — Meteoroid Fluxes at the Formation of OP f MT f 
and LH 


Area 

Age (10» yr) 

F (m, t) 

OP 

3.3 

5.3X10 3 4> 

MT 

3.65 

1.3X10 4 

LH 

4.6 

1.6X10 5 


according to Shoemaker (Gault, 1970). Gault's 
conclusion, assuming a time-constant influx, was 
consistent with the crater distributions in MT. 
However, a time-constant influx cannot explain 
crater distribution in the Lunar Highlands or 
explain the difference in distributions between 
MT and OP related to their absolute formation 
ages. But these can be explained if one accepts 
the flux to be 

$=10 - 14 *mr l [m- 2 s- x ] (7) 

This flux function, equation (7), is lower than 
the one given by Naumann, Hawkins, and Shoe- 
maker. But equation (7) does agree with the flux 
given by Watson (1956) and will be discussed in 
the next section of this paper with respect to 
other flux data. 

If one accepts, for the moment, this flux func- 
tion to be valid for today, the time-variable flux is 

F(m, t) = 10 -I4 - 6 m _I e -2 * 6 * (8) 

The comparison with the crater-number densi- 
ties in MT and OP works if one further accepts 
the D/d values for large craters to be on the order 
of 10. This value is quite low, but the correspond- 
ing value of 2 for /x-sized craters is experimentally 
known (Bloch et al., 1971; Mandeville and 
Vedder, 1971). This information suggests that 
the flux assumed is still too high for the large- 
crater range (km diameters). The fact that the 
meteoroids show a d~ z distribution and the craters 
show a D -2 ' 9 distribution is interpreted in the 
variability of the D/d values of 2 for jt-sized 
craters at 20 km/s impact velocity and increasing 
with increasing meteoroid diameters. 

The distribution of smaller craters ( < 140 m 
diameter in MT and <50 m diameter in OP) in 
the diagram of figure 6 shows a D ~ 2 distribution 
for MT and OP. As interpreted by Shoemaker 
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et al. (1970) and experimentally shown by Gault 
(1970), small craters are lacking because of 
meteoroid impacts (for example, from erosive 
effects of ejecta material), solar wind (from 
sputtering effects), and possibly thermal gradient 
influences. This distribution could be verified 
down to crater diameters of about 50 a as shown 
by Horz, Hartung, and Gault (1971) and by 
Neukum (1971). 

The time-variable formalism can be applied to 
crater sizes in the equilibrium status and mean 
crater lifetimes t can be calculated by the formula 

N 

e~ 2M dt= - 

where 

N/F= cumulative number of craters on the 
equilibrium curve 

$(m) — cumulative flux of the corresponding 
meteoroids. 

For 0.3 mm diameter craters (observable on 
small lunar surface samples) the mean crater 
lifetime is approximately 10 7 yr. For 1-m diameter 
craters the lifetime is 10 8 to 10 9 yr which is con- 
sistent with cosmic-ray-exposure age determina- 
tions for the Apollo 12 soil of 350 million yr 
(Kirsten, Steinbrunn, and Zahringer, 1971). 

The latest progress has been gained by scan- 
ning smooth-crater glass linings of sample No. 
12063.106, glass-coated sample No. 12024.8.1, and 
a glass ellipsoid from sample No. 14257 (Feehtig, 
Mehl, Neukum, and Schneider, 1972). It was 
possible to discover and count craters down to 
about 0.3m diameters with a D -08 size distribution. 
The crater number density for craters >1m 
diameter is found to be 10 3 craters per cm 2 . The 
conversion into fluxes according to Neukum and 
Dietzel (1971) results in a m~° 6 distribution down 
to 0.15m diameter particles (fig. 7) . 

The overall meteoroid flux distribution as 
derived from the distribution of lunar craters is 
shown in the flux diagram of figure 8. The general 
slope in the mass distribution is — 1 according to 
equation (7). The distribution shows a depiction 
at about 30M-diameter particle size. There, the 
slope changes to -0.6 down to meteoroid masses 
of about 10 -15 g. These submicron-sized craters 
are considered as an important observation con- 
cerning the existence of submicron-sized solid 


particles in the interplanetary space. The lower 
slope in the mass distribution may be interpreted 
by a variable cutoff for small particles due to the 
Poynting-Robertson effect. This variability is 
caused mainly by different eccentricities of the 
particle orbits (Dohnanyi, 1969) and a variability 
of the chemical composition of the particles. 
According to Shapiro et al. (1966) metallic 
particles show a much higher cross section for the 
influence of the solar electromagnetic interaction 
than nonconductive particles. Therefore, the 
smallest particles would be the nonmetallic 
particles since the metallic component cannot 
survive as close as 1 AU to the Sun. 

THE PRESENT METEOROID FLUX 

Figure 9 shows the flux of cosmic dust cumula- 
tively plotted as a function of mass. The results 
are discussed in four groups: 

(a) Results from lunar crater distributions 

(b) Results obtained in interplanetary deep 
space (including zodiacal light measurements) 

(c) Spacecraft results (between Earth and 
Moon and from Lunar Orbiter) 

(d) Results from dust collectors and detectors 
in the atmosphere (rocket- and balloon-borne 
experiments) . 

Lunar Craters 

The general mass distribution derived from 
lunar-crater statistics have been discussed in the 
previous section. The results are labeled “MOON” 
in figure 9. The slope for larger particles (>30m 
diameter) is —1, for smaller particles (<30m 
diameter), —0.6. 

Interplanetary Deep Space 

Measurements were reported by Alexander et al. 
(1971) from the Mariner IV dust experiment and 
by Berg and Gerloff (1971) from the Pioneer 8 
and 9 experiments. Leinert (1971) has calculated 
fluxes from earlier zodiacal-light measurements by 
Elsasser (1958), Ingham (1961) and Weinberg 
(1964). 

From the Pioneer dust experiment a cutoff was 
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Figure 7. — Microcrater distribution on various lunar samples. 
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Figure 8. — Cumulative meteoroid fluxes as a function of mass. 


reported (Berg and Gerloff, 1971 j Gcrloff and 
Berg, 1971) at about 10“ il g and interpreted with 
orbit elements according to Dohnanyi (1969). 
Generally, the in situ results and the indirect 
measurements of the zodiacal light are very 
similar. The cutoff from the zodiacal-light meas- 
urements is a question of interpretation. 

A comparison of these results with the results 
from lunar-crater statistics leads one to question 
the Pioneer detection cutoff of particles at 10 -11 g. 
On the contrary, sub micron-sized lunar craters 
have been found, thus indicating the existence of 


submicron-sized particles in the interplanetary 
space. Two possible explanations of this apparent 
discrepancy are as follows: 

(a) If a variable particle-size cutoff exists 
because of chemical composition, the smallest 
particles are nonmet allies (silicates, quartz) . The 
general lack of metals for the smallest particles 
influences heavily the ion production of the impact 
detectors used in the Pioneer mission as discussed 
in a previous section of this paper. 

(b) The Pioneer dust experiment with a quasi- 
circular orbit around the Sun cannot detect 
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DUST METEORS 



Figure 9. — Cumulative cosmic dust fluxes as a function of mass. 


circularly orbiting particles, since the relative 
velocity is not high enough for small particles to 
penetrate the thin foil in front of the detector and 
produce enough ions to be recorded. (The lowest 
relative velocity must exceed 0.7 km/s to be 
detected.) 

Spacecraft Results 

Higher fluxes have been measured by Earth 
and Moon satellites. Using different techniques, 


results were reported from the following experi- 
ments: capacitive detectors on Pegasus (Anon., 
1966, 1967), penetration-experiment on Explorer 
XXIII (Naumann, 1966), microphones on Cosmos 
213 (Nazarova and Rybakov, 1971), ionization 
detector on OGO III (Alexander, Arthur, and 
Bohn, 1971), ionization detector on Lunar 
Orbiter 35 (Alexander et al., 1971), and the 
S 10- S 12-craters (Hemenway et al., 1968). Two 
possible explanations for these higher near-Earth 
fluxes are : 
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(a) Temporarily higher influxes of lunar debris 
are produced both by impacts of meteor streams 
and sporadic dust particles on the lunar surface. 
This origin was discussed by Alexander et al. 
(1971). Gault et al. (1963) and Colombo et ah 
(1966) have also discussed this problem. 

(b) Fireballs that might be of low-density 
material (Dubin, 1971) form disintegrated prod- 
ucts. 

Atmospheric Dust 

The highest flux numbers are reported from 
experiments in the atmosphere using rocket and 
balloon-borne dust collectors and detectors. Based 
on the Pandora and the balloon-top collections by 
Hemenway et al. (1971), the balloon collections 
of Brownlee et al. (1971), the MPI 68, 69, 70 
collections (Auer et al., 1970; Fechtig and Feuer- 
stein, 1970; Fechtig, Feuerstein, and Rauser, 
1971), the LUSTER 68 flight through a NLC 
display (Farlow et ah, 1970), and the detector 
results of “Lund” (Lindblad et ah, 1970), MPI 
70, 71 (Fechtig, Feuerstein, and Rauser, 1971; 
Rauser and Fechtig, 1972), one can state the 
following items : 

(a) No cutoff down to masses of about 10 -15 g 
could be found. 

(b) Easily fragile particles are evident and 
support Jacchia’s (1955) and Verniani’s (1969) 
observations of low-density material. 

(c) Large variations probably due to meteor 
show T cr activities and Noctilucent Cloud displaj^s 
are observed for the fluxes. 

(d) Still not sufficiently known is the influence 
of the upper atmosphere on the dynamics of small 
particles. In 1970 w r e (Rauser and Fechtig, 1972) 
measured a velocity profile of micron-sized 
particles between 70 and 110 km altitude. 

The measured velocities are low r er than expected 
at 110 km altitude which suggests that the 
particles arc of a fluffy low-density material. 
However, the measured velocities are about 2 


orders of magnitudes higher than expected (Korn- 
blum, 1969) at 70 to 80 km altitudes. The exist- 
ence of a temperature minimum at about 85 km 
altitude leads one to expect a diameter and/or 
density increase of particles by absorption of 
condensable material like w r ater or C0 2 . This 
mechanism might also explain the Noctilucent 
Clouds. A layering effect can be calculated from 
the detailed data. 

Finally, total influx rates can be calculated from 
the flux curve. One calculates a total daily influx 
of interplanetary material on the Earth of about 
25 tons. This number is consistent with investiga- 
tions by Keays et al. (1970). These authors have 
investigated lunar samples for meteoritic com- 
ponents. Although the investigations were carried 
out in a completely different scientific field, the 
results show an agreement. 

SUMMARY 

This paper contains a description of techniques 
used in recent experiments to detect and analyze 
cosmic dust and micrometeorites. Furthermore, 
the results both from the study of lunar crater 
statistics and from in situ measurements have 
been reviewed. 

Ionization detectors and controlled collectors 
represent important progress in the techniques 
used in the research field of cosmic dust. 

The results from lunar crater statistics show an 
agreement with the results obtained from in situ 
measurements in interplanetary space and derived 
from zodiacal-light measurements. The near- 
Earth results show an enhancement in the flux 
numbers. This can be caused either by secondary 
lunar debris or by disintegration of low r -density 
fireballs in the outer atmosphere. 

For future experiments in situ measurements are 
needed in interplanetary space. The question of 
the chemical composition of particles is important 
with respect to the origin of cosmic dust. 
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